Numerical Evaluation of Shot Noise using Real Time Simulations 
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We present a method to determine the shot noise in quantum systems from knowledge of their 
time evolution - the latter being obtained using numerical simulation techniques. While our ultimate 
goal is the study of interacting systems, the main issues for the numerical determination of the noise 
do not depend on the interactions. To discuss them, we concentrate on the single resonant level 
model, which consists in a single impurity attached to non-interacting leads, with spinless fermions. 
We use exact diagonalisations (ED) to obtain time evolution, and are able to use known analytic 
results as benchmarks. We obtain a complete characterization of finite size effects at zero frequency, 
where we find that the finite size corrections scale oc , G the differential conductance. We also 
discuss finite frequency noise, as well as the effects of damping in the leads. 

PACS numbers: 



The study of current fluctuations in nanodevices such 
as quantum point contacts and tunnel junctions is deeply 
connected with some of the most important physical 
questions. These include the nature of fundamental ex- 
citations in strongly interacting electronic systems [l|, Q , 
the possibility of fluctuation theorems out of equilibrium 
3|, and the time evolution of many body entanglement 
4 , m . Experimental progress in this area has been swift 
- second and third cumulants have been measured in sev- 
eral systems 0,0] , shot noise of single hydrogen molecules 
have been measured [1] , and even the full counting statis- 
tics has been obtained in semi conductor quantum dots 



On the theory side, the free case has given rise to a 
lot of analytical studies d, [l^, EH, but progress on 
the most interesting situations - far from equilibrium and 
with strong interactions - has been difficult (see [l^l for a 
review). Over the years, extensions of the Bethe ansatz 
to study transport properties have been proposed [isl - 
[l6j . which might open the road to important progress. 

The numerical situation is somewhat similar. I — V 
characteristics had remained inaccessible until very re- 
cently, where methods based on time dependent DMRG 
have been successfully used to tackle some simple nanos- 
tructures [l7| . 

It seems reasonable to expect that time dependent 
DMRG results could also be used to determine current 
fluctuations, which could also, in some setups, be de- 
termined analytically [l8| . In order to reach this ambi- 
tious goal, it is crucial to be able to extract cumulants - 
in particular the shot noise - from real time simulation 
methods. We propose in this paper a method to do this. 



The main problem in the determination of the noise 
- and the main emphasis of our work - is the finite size 
analysis of the results of non-equibrium correlation func- 
tions for finite systems. To concentrate on this aspect, we 
only discuss results for the non-interacting resonant level 
model (RLM) where the numerical data can be obtained 
using exact diagonalisation (ED) techniques. Since in 
this specific case there are straightforward analytical so- 
lutions of the problem, we can check in great detail the re- 
liability of our approach. We want to emphasize however 
that the method to be described below is independent 
of the underlying numerical algorithm, and can therefore 
also be used in combination with, e.g., time dependent 
Density Matrix Rcnormalisation Group (td-DMRG). We 
will report on this in the case of the interacting resonant 
level model in a forthcoming paper. 

Wc note that prior to our work, a numerical study 
of the full counting statistics for another non interact- 
ing model appeared in [l^. The method used there is 
however entirely specific to the free case, and uses inter- 




FIG. 1: Sketch of the IRLM model. J' denotes the coupling 
of the impurity to the left and right lead, J is the hopping 
element in the leads, U is the interaction on the contact link, 
and is a gate voltage. In this work we set = (the dot 
level is on resonance), and (7 = (non- interacting case, which 
allows to treat the problem via ED). 
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mediate analytical results from [l^l- Our approach, in 
contrast, is based directly on the 'experimentally' mea- 
sured time evolution of the current. 

To make things concrete, we start by giving the Hamil- 
tonian of our test system. It is composed of a structure 
(i?s) attached to leads (Hl) described in real space by 

H = Hi^ + Hs, (1) 

— 2 oo 

= ~J ^L+l^m ~ ^ ^L+l^m + h.C, (2) 

— oo m— 1 



Hs = -J' ("^'"^0 + 4cm) + Vgno 



Tn=±l 



m=±l 



(3) 



cf. Fig. [TJ The nearest-neighbour hopping matrix ele- 
ments in the leads and the coupling of the structure to 
the leads are given by J and J'. In the remainder of 
this work we concentrate on the resonant case at zero 
gate voltage Vg = and half filling. Since we want to 
compare the numerical data with analytical results, we 
furthermore restrict ourselves to the non-interacting case 
with [/ = 0. 

To prepare the system in a state with finite current 
through the structure, we add a charge imbalance opera- 
tor Q = Vsu (^Ni^ — Njij /2 to the Hamiltonian and calcu- 
late the initial state as the ground state \'^{t = 0)) = |^'o) 
oi H + Q. Here, TVl (-/Vr) counts the particle number in 
the left (right) lead. We then perform the time evolution 
with the Hamiltonian without Q. The time evolution 
is performed by means of the time evolution operator 
U = c~'^*, while all expectation values are evaluated 
with respect to the initial state |*I'o)- For details see [20l - 

m. 

The current operator /,„ for the current at bond m is 
given by 

In (t) = Jm [c] (OS + l (t) - c]+i (t)] ■ (4) 

We define the current operator as an average over the 
current on the left and the right contact of the nanos- 
tructure 
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FIG. 2: Time dependent current I{t), and current correla- 
tion function S{t, tmin) with tmin ~ 12, iu the non-interacting 
resonant level model RLM, with tight-binding leads and a 
finite system size of M — 60 lattice sites, for different val- 
ues of the bias voltage Vsd- The I{t) curves show the three 
time regimes given by the settling time ts and the recurrence 
time t^. The highlighted time domain indicates the integra- 
tion range [t 



mill , "^max 



The inset demonstrates an additional 
subtility: the correlation function shows finite size reflection 
effects on the time scale t — tmin ^ iR/2, which imposes an 
additional restriction on tmax- 



current correlations in the time domain 

S{t,t') = i(A/(t)A/(i') + A/(t')A/(t)) (6) 



Re(A/(t)A/(t')> 



(7) 



have to be calculated in a non-equilibrium zero- 
temperature state, where A/(t) = /(t) - (/(t)) [ii,!!^]. 
Therefore, the time dependent expectation value 

(A/(t)A/(t')> = (4'o|c'^*A/c-'^(*-*')A/e-'^*'|*o) 

(8) 

has to be evaluated. In a steady state the correlation 
function must fulfil S{t,t') = S{t — t'). Then the noise 
power can be defined as the Fourier transform 

2'n:S{uj + uj')S{uj) = (A/(w)A/(w') + A/(w')A/(a;)), (9) 

where 

oo 

S'(a;)==2 / die'"*S'(t,i' ==0) (10) 



(5) 



The expectation value of /(t) in the RLM for J' = 0.4J 
and for some values of 14d is shown in the upper part 
of Fig. [21 Effects like the finite settling time ts and the 
finite transit time tR. as well as the I-VJjD-characteristics 
have been discussed before in great detail [1^ [13, HI] . 

Shot noise is defined as the zero-temperature contri- 
bution to noise in a transport state. To obtain the noise 
power spectrum from a real time simulation, the current- 



= 4Re J d<e'"*5(i,t' = 0). (11) 



The right-hand side of the equation accounts for the sym- 
metry S{t — t') = S{t' — t). In a steady state, of course, 
this expression should be independent of the choice of the 
time t' 



oo 

S = 4Re j dte"^^*-^'^S{t,t') Vt'. 



(12) 
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FIG. 3: Noise S and squared differential conductance for 
J' = 0.4J. The blue lines represent the analytic values ob- 
tained using the Landauer-Biittiker approach. The finite size 
of the system introduces an additional noise proportional to 
G^/M (for discussion see text). 
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FIG. 4: Noise S{uj) vs. frequency uj, both rescaled with re- 
spect to the width of the conductance peak F = (2J' /J)^ , for 
different values of the bias voltage. The lines going through 
the numerical values (represented by crosses) are just guides 
for the eye. The dotted lines correspond to the analytical 
result for the wide band limit. 



In the zero-frequency limit a; = this expression simpli- 
fies to 



S = S{uj^0)^4: / dtS{t,t') 



oo 

4 J dtRe{Ai{t)Ai{t')). 



(13) 



(14) 



We now want to see whether the noise can be reliably 
obtained in a real time simulation based on this formula. 

There are of course many obstacles. The first comes 
from the calculation of a nonequilibrium correlation func- 
tion in the time domain from a real time simulation. Be- 
cause we are restricted to a finite system with AI lat- 
tice sites and hard walls, a steady transport state is not 
well defined. Instead, we make the attempt to calcu- 
late the time evolution from the initial non-equilibrium 
state \'^{t = 0)) as described before and look for a 
quasi-stationary time regime. The "switching" of a finite 
source-drain voltage Vso at initial time causes a ringing 
of the current [ll], cf. also Fig. [21 which decays expo- 
nentially within a settling time ts- The current finally 
enters a plateau regime, where the size of the plateau is 
given by the recurrence time which is finite due to the 
finite size of the system [2l[, Fig. [2l 

Having obtained results for the system "close" to a 
steady state, to obtain the quantity S we now have to 
evaluate the integral (fT4|) in a limited time range 



dtRe(A/(t)A/(t,„i„))* 



(15) 



where t^i^ > to, and tmax < ^r- In a hypothetical sit- 
uation with a system of infinite size where tR — oo the 
contribution of dtRe{AI{t)AI{trnin))^ can be ne- 
glected if Re(A/(t)A/(tmin))* is small for t > tmax as 



compared to the mean value in the range imin < t < tmax- 
One therefore has to choose the size of the system big 
enough to ensure the correlation function to drop to zero 
within the recurrence time. 

The finite recurrence time introduces a finite cut- 
off frequency Wcut = Stt/^r cx 1/M. This is the main 
problem we encounter. In contrast to the situation of 
infinite leads, where zero frequency noise vanishes with- 
out applied voltage, we now find a contribution to the 
zero voltage shot noise of the order of ^(wcut)! The low 
frequency domain is the most interesting for the kind of 
problems we wish to study: low frequency is low energy 
and thus strong coupling between impurity and leads. 

The magnitude of the finite size effects for the type of 
systems that can be studied today is far from negligible. 
On fig. [3] we give results for the shot noise ^num obtained 
for different system sizes of M lattice sites, as well as the 
expected result S in the thermodynamic limit obtained 
from the Landauer-Biittiker approach (this is discussed 
in more details in the appendix). ^^W?'^^//^/ While the 
results measured for finite size and the asymptotic results 
agree at large voltages, there is a marked difference at 
small voltages, with an offset at vanishing Vsd- On the 
figure, we also represented the finite size correction 



A5„ 



s 



(16) 



rescaled by the system size M. For different values of M 
the rescaled finite size corrections M x A^num collapse 
very well on a single curve, indicating that the main fi- 
nite size effects scale linearly with 1 /M in the considered 
parameter regime. One may expect that the cut off given 
by the finite size of the leads corresponds to an effective 
finite temperature ~ resulting in a low voltage off- 
set ^ G/M. However, we find AS'num cx; with the 
differential conductance G{Vsd) = dl (Vsd) / dVso ■ 

To understand the behavior of AS'num, we consider the 
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FIG. 5: Slope of the frequency dependent shot noise in the 
limit tj — >■ 0^, rescaled to fit with G^, for different system 
sizes. In the low voltage regime we find finite size effects. 

full frequency dependence of the shot noise. It can easily 
be obtained analytically in the wide band limit - see the 
appendix, Eq. (jA6p . For values of J' / J ^ 1, the numer- 
ical results obtained for the model with cosine dispersion 
relation should be consistent with the analytical result 
as long as the considered frequency is small compared to 
the band width. This is illustrated on fig. ID There, the 
frequency dependent noise is obtained via 

^„um(c^) = 4Re j dte'"(*-*-")5(t,W) (17) 

tmin 

for different values of the bias voltage Vsd . For big values 
of w, the effects of the band curvature arc quite marked 
- as can be seen by the departure of the various guide 
lines from the dotted lines representing the analytic wide 
band limit results . 

To understand the voltage dependency of the finite size 
corrections, we consider the low frequency behaviour of 
the analytical results in the wide band limit where we 
find 

S{uj > 0, Vsu) = ^(0, Vsu) + A5(w, Vsu) + Oiuj^) (18) 

with the correction in first order with respect to uj 

A5(w,Fsd)ocG2(Fsb)l^. (19) 

For the system with finite band width, we have checked 
this expression by extracting the slope dS{uj, Vsu)/du} in 
the limit cj — > 0^ from the numerical data. Again we 
find good agreement with in a voltage regime where 
finite size efi^ects can be neglected. Fig. O 

Inserting the cutoff frequency now leads to the expres- 
sion 

A5(weut,^sD)« ^G2(y,„) (20) 
which is in good agreement with AS'num(14D), cf. Fig. [31 



FIG. 6: Noise S and squared differential conductance 
for J' = 0.4J. The blue lines represent the analytic values 
obtained using the Landauer-Biittiker approach. The system 
size is fixed to M = 60 lattice sites, while at the boundaries, 
the hopping matrix elements are exponentially damped with 
the damping constant A"^''^ on Ma links. This results in an 
effectively enlarged system with Mca lattice sites. The finite 
size correction ASnum, here rescaled by M^s, again collapses 
on a single curve for different MeS, and is still proportional 
to G^. 

Using our knowledge of the finite size correction, we 
can now control the extrapolation of numerical data: in 
Fig. [7] we show the results obtained using linear extrap- 
olation 1/M ^ for J' = 0.3 J and J' = 0.4 J. We find 
indeed very good agreement with the analytical result. 

The non-interacting case is of course very simple to 
calculate numerically (regardless of the possibility of the 
Landauer-Biittiker treatment). The numerical main ef- 
fort consists in the exact diagonalisation of the M x M 
Hamiltonian matrices as well as the calculation of the 
time evolution which involves the multiplication of M x 
M matrices. Including interaction spoils this approach. 
Instead, one has to resort to approximative time evo- 
lution schemes using methods for correlated electrons. 
In [i^l we showed that using Wilson leads, or damped 
boundary conditions (DBC), respectively, with a weak 
damping constant allows one to effectively increase the 
system size to M^s > M lattice sites without changing 
A'/, where a rough estimate for Afofr has been given as a 
function of the damping constant A and the length of the 
damped leads Ma 

Meff « M - 2A/a + ^ (A^^^/2 _ (21) 

We now use this estimate to perform the linear extrapola- 
tion to infinite system size, where we additionally adjust 
the estimate by fixing the extrapolated value to analytic 
results (cf., e.g., [3]| ) 

S{Vsu = 0) = 0. (22) 

To verify this approach we performed calculations for a 
non-interacting system with M = 60 lattice sites and 
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FIG. 7: Shot Noise S as function of the source drain vohage 
VsD in the non-interacting resonant level model. The analyt- 
ical result was obtained using the Landauer-Biittiker theory, 
Eq. (|A5[I . while the numerical result is computed for systems 
of different finite sizes of M = 120. ..180 lattice sites with a 
subsequent linear extrapolation of 1/M — >■ 0. The two curves 
correspond to different couplings J' of the impurity to the 
leads. Furthermore, we used DBG in order to effectively in- 
crease the system size. Here, the system size was fixed to 
M = 60 lattice sites. For weak damping and for not too big 
values of Vsd, we find very good agreement with the undamped 
case. 



DBC, for J' = 0.4 J. For the damped leads we used differ- 
ent combinations of A and Ma, where we used values for 
the damping constant in the range g [0.93, 1.0] for 

damped leads of A/a = ... 26 lattice sites (while keeping 
the total number of lattice sites M fixed!). The estimate 
for the effective system size, Eq. pip, is checked by look- 
ing at the scaling behaviour of the finite size correction 
AS'num; where we now find linear scaling cx 1/Afcff, cf. 
Fig. H 

The result is shown in Fig. [71 We find remarkably 
good agreement with the analytical result, while we have 
to point out that, for values of the bias voltage in the 
order of the band width, the approach fails, which has 
to be expected since the estimate of the system size only 
works in a limited voltage range (27| . Additionally we 
find the numerical data to be very noisy depending on 
the respective configuration of the damping conditions. 

To conclude, we have introduced a new way of extract- 
ing the finite bias shot noise from real time evolution cal- 
culations. Very accurate quantitative results are possible 
as long as finite size effects arc treated properly. The 



possibility of effectively enlarging the system size using 
Wilson leads has been successfully tested for a limited 
voltage range. Our results for shot noise in the RLM 
show very good agreement with results obtained from 
the Landauer-Biittiker approach for zero as well as fi- 
nite frequencies. Additionally we confirm a depen- 
dent scaling of the finite size error in the zero frcqunccy 
regime. 

The concept is not restricted to non-interacting 
fcrmions and could be implemented using numerical 
methods for interacting quantum systems. In this case 
as well, we expect the finite size corrections to go as 
1 jM because the cutoff frequency Wcut has the same de- 
pendence. We note however that the prefactor might 
not be G^(Vsd) exactly: the jury is still out on the 
small frequency behavior of the noise in the presence of 
interactions, [s^, HI] 

Wc acknowledge the support by the DFG Center for 
Functional Nanostructures (CFN), project B2.10. 



Appendix A: Analytical results 

For the non-interacting system in the thermodynamic 
limit, Eqns. (IT][3]), with [/ = 0, one can derive analytic 
results for the differential conductance and the shot noise 
from single particle scattering states. The energy depen- 
dent single particle tunneling probability r is given as 



1-6V(4J2) 



1 + (e ~ Fg) ( J2 - 2J'2)e - jVg /(4J'4) 



(Al) 



Using the Landauer approach (cf., e.g., [31j ) one obtains 
for the dc current / and the noise S at zero frequency 
and for a symmetric voltage drop Isd 



/("K,d) 
5(14d) 



Vhd/2 

-v:,d/2 



der(e). 



der(e)(l-r(e)). 



For the resonant case with = this results in 



G(ysD) 



1 - vJkx^P 



dl 

dVsr, " 1 + Vsu\j^ - 2 J'2)/(16J'4) - 



(A2) 
(A3) 

(A4) 
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where Tb = AJ'^ /\/ ,P — 23'"^ is a scale. For the frequency dependent noise, we will contend ourselves by giving the 
result for the wide band limit only: 



S{i^, Fsb,Tb) =^e(FsD - \uj\)\^ [arctan (-|^) + arctan 



^sD - 2|w| 



I arctan ^ 



V^SD + 2|^| 



arctan 



Tb / T| + (ysD-2|c.| 



-)} 



(A6) 



4 / 

where in this limit — -y- 
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